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ON REGULARITY OF THE TIME DERIVATIVE FOR DEGENERATE 

PARABOLIC SYSTEMS 

JENS FREESE AND SEBASTIAN SCHWARZACHER 


Abstract. We prove regularity estimates for time derivatives of a large class of nonlinear 
parabolic partial differential systems. This includes the instationary (symmetric) p-Laplace 
system and models for non Newtonien fluids of powerlaw or Carreau type. By the use of 
special weak different quotients, adapted to the variational structure we bound fractional 
derivatives of Ut in time and space direction. 

Although the estimates presented here are valid under very general assumptions they are 
a novelty even for the parabolic p-Laplace equation. 


1. Introduction 

In this paper we will prove fractional differentiability for the time derivative ut- This we 
can do for a very general class of essential non-linear PDE; the model case is the parabolic p- 
Laplace. But our results include models of non-Newtonian fluids of power-law type initiated 
by Ladyzhenskaya [20l |2T] and J.J.-Lions [23] (see Subsection 11.31 below). 

To our knowledge little regularity for the time derivative is known up to now even for the 
homogeneous instationary p-Laplace equation (he. dni), with / = 0). In comparison quite a 
lot of regularity is known for the space gradient. After the pioneering work of Friedman and 
DiBenedetto |7| where they proof that space gradients are Holder continous if the right hand 
side is zero a multifarious collection of results on gradient regularity was developed. We will 
discuss the known results and its implication of the regularity of ut more detailed in the next 
subsection. 

In the framework of embedding theory, interpolation to finite element spaces or com¬ 
pactness issues the regularity of the time derivative is of essential importance. Indeed, its 
regularity often restricts the estimates crucially. Therefore and because of the generality of 
our approach, we emphasize that the results and techniques developed in this article will be 
an important step in the analysis and numeric to many instationary applications. 

1.1. The parabolic p-Laplace system. The reader will be introduced to the results of this 
article by explaining them on the model case. Namely, the (symmetric) parabolic p-Laplace 
system. 

Ut — div{\Duf~‘^Du) = / on Qt 

(1.1) u = 0 on (0,T) X dQ 

u{0,x) = uo{x) on n. 

Here H C M"’ is a domain in space and [0, T) a time interval Qt ■= [0, T) x Q. We will 
use Dv as a substitute for both the gradient Vu as well as the symmetric gradient ev := 
^(Vu -|- (Vu)"^ (which is of course only defined, whenever n = N). This is due to the fact, 
that we can treat both cases simultaneously. Standard existence theory implies that there 
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is a unique solution u € LP{{0,T),Vq’^{^)) n L°°{[0,T), L‘^{Q)) for every right hand side 

/ € LP'{[0,T), and every uq G L‘^{n). The space := . 

Our aim are estimates for ut in terms of Nikolskij spaces. These are spaces that represent 
fractional derivatives. For a G (0,1] (the order of derivative) and q G [l,oo) (the exponent 
of integrability) we say that g G AA"’'?((a, 6 )), if 




,H{a,b)) 


:= sup 

/i>0 


rh—h 

J a 


g(t + h)- g{t) 


h° 


dt+ \g\'^ dt < oo. 


For a general domain A we say that g G if the property above holds in all coor¬ 

dinate directions. Nikolskij spaces are closely related to fractional Sobolev spaces VF“’'^(^). 
Let us just mention that 0 < a < /3 < 1 both 

lF^’''((a,6)) C A^^’‘'((a,6)) C W^’'^{{a,b)). 


In our setting the relation can be used, as for fractional Sobolev spaces the following embed¬ 
ding theorem is available: For a Lipschitz domain A C M” and g G W°‘’‘^{A) we have 


( 1 - 2 ) 

This implies for instance that if <7 G J\f°‘’'^{A), g G for every 1 < a < For the 

embedding theorem and a more detailed study on the given function spaces we refer to [291 [2]. 
For the study of fractional spaces in the framework of PDF we recommend [25]. 

The hrst result of this article is that ut has fractional derivatives of order one half in time. 
More precisely ut is in the Bochner-Nikolskij space T — a), L^(n)); which means 

sup / / —^dxdt < oo, 

h<a . 


r° f 

Ut{t + h) - Utit) 

la J Q. 

hh 


for a > 0 (see Theorem ESI 

Our second result concerns spatial derivatives of ut- In the degenerate case the existence 
of the space gradient of ut is not necessarily available. Therefore the fractional regularity of 
Ut in space direction is of interest. 

In case of the whole space hi = M"' (or in the space periodic case) we can prove that ut has 
fractional spatial derivatives. 

For all 1 < p < oo we prove a global estimates, i.e. 


sup 

h<a 


r° / 

Ut{t,x + het) - Ut{t,x) 

la Jmp 

hi 


2 

dxdt < oo. 


for Ci the i-th unit vector. This implies that ut G L^((a, T —a), AA4’^(M”)) (see Theorem l2.4ll . 

These estimates can be refined in case the solutions are space-periodic. In this case we 

can show that ut G (I^”))) provided 2 — ;| < p < 4, furthermore we obtain 

Nikolskij differentiability orders in the interval ( 5 , 5 ] for all 2 — ^ < p < 00 (see Theorem l2.6l) . 

It turned out to be more difficult to gain regularity in space direction then in time. If one 
considers bounded domains one gets non-trivial additional boundary terms; a possible way 
of treating also this situation looks sophisticated and we postpone it to a future project. 

In the special case of the parabolic p-Laplace one will find quite satisfactory estimates of 
u and Vrt in the literature. Particular, the gradient is Holder continuous provided the right 
hand side is smooth [Gj Chapter 8 ] and [261 [2^ . This is the parabolic version of the elliptic 
analogue initiated by Uhralzeva m for equations and Uhlenbeck [30] for systems. Moreover, 
an L'^-Theory [T] (by non-linear Calderon-Zygmund theory initiated by Iwaniec [15] 1 and 
pointwise estimates of Riesz type are available for the parabolic p-Laplace [T9l[T8|. The only 
results for the time derivative that is available is see [22] and [231 Theorem 8.1]. 

In the next two paragraphs we introduce two generalizations of the parabolic p-Laplace 
for which we can prove the very same regularity results. 
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1.2. General systems. We present the general case of systems where the elliptic operator 
is Euler of a variational problem. 

We consider the following energy functional. For u : —>• 

(1.3) J{v)= / F{x, Dv) — fvdx. 

Jn 

The respective Euler operator is defined as —div(5^E). Again we will use Dv as a substitute for 
both the gradient as well as the symmetric gradient as we can treat both cases simultaneously. 

At this point we discuss the assumptions for systems with p-growth. Observe, that we will 
include later generalized conditions of Orlicz growth (see Assumptions [22] below). 

Assumption 1.1. For x G M" and Q,P & we assume the following growth conditions 
for ^ > 0 

(a) F is measurable and F{x, ■) G C'°’^(]R”'^^) 0 \ {0}). 

(b) For Aij{x,Q) := ^Q^^F{x,Q), A(|(5|^ + h'^)^\Q\'^ < ■ Q and \A{x,Q)\ < 

A(|(5|^ + ff^)^\Q\, for 0 < X < A. 

(c) {A{x, Q) - A(x, P)) • (g - P) ~ \V{Q) - V{P)\\ for V{Q) := (^^ + |g|2)^g. 

(d) \D^A{x,Q)\<{\Q\^ + g.‘^Y-^\Q\. 

(e) DqA{x, Q)P®P^ |P|2(/i2 + 

(f) F{x,0) =0. 

In case we consider the symmetric gradient, we need additionally that 

(g) If Q is symmetric, then A{x, Q) is also symmetric. 

The constants used in this assumption are called the characteristics of F. 

Remark 1.2. Well known is the case of Uhlenbeck structure; F{x, \Du\). In case F{x, Du) = 
^\Du\^ (ll.ip becomes the p-Laplace or the symmetric p-Laplace. 

It is possible to replace (b) by the assumption 

A(|g|^ + — k{^) < Q) • Q Q)\ < + k{x), with 

k,K e L^{W^)r\L°°{W^). 

It is also possible to replace (d) by the assumption DxA{x,Q) = Bi{x,Q) + B{x), such 
that |ili(x,g)| < c(|g|^ + fi‘^)^\Q\ and VB G Lp'{Qt). More general one considers B to 
be of lower order. I.e. of the form B{x, Q) has lower then p — 1 growth in \Q\. Finally it is 
possible to replace (f) by F(x,0) = g{x) G L^{Qt)I^L°^{Qt)- bhe did not include this weaker 
assumptions, as the proof would not change significantly but would be more complicated to 
read. 

From the above energy functional (jl.3p we associate the following flow. 

ut — div(A(x, Du)) = / on Qt 

(1.4) M = 0 on (0,r) X dH 

u{0,x) = uo{x) on n. 

In the general case, in particular problems where the symmetric gradient is involved one knows 
only the improved differentiability based on classical difference quotient techniques. Recently 
in the PhD thesis of Burczak, some local (in space) estimates of fractional time derivative of 
u is shown [5] Lemma 4.9]. See also [1] for an overview of regularity of symmetric p-Laplace. 
Our paper gives an additional progress in the theory, by showing partial differentiability of ut 
in time and space direction. Clearly there are many paper on partial regularity or short-time 
regularity, which is not of our concern here. 
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1.3. Non-Newtonian fluids. Our regularity theory implies results for certain non Newto¬ 
nian fluids; in explicit for solutions of the following system 

ut — div( 74 (x, eu)) + Vtt = / on Qt 
div(tt) = 0 

fl ^ 

n = 0 on (0, T) x dO. 

n( 0 , x) = no(x) on fl. 


Here tt : 0 —)• M” describes the velocity and vr : w —>■ M describes the pressure of a fluid 
moving through time in fl. We denote the symmetric gradient as ev := ^(Vu -|- (Vu)^). The 
vector field / : 0 —)• is the outer force (e.g. the gravity) and uq some given initial state. 

The matrix A{x, eu) is the (non-linear) stress tensor. Our model situation is, that A is of 
power law type, i.e. A{x, eu) = -|- \su\)^^eu, for 1 < p < oo, v > Q and ^ > 0. However 

we can also treat Carreau type fluids: A[x,eu) = Voo^u + -|- \eu\)~eu, for Uoo > 0. 

Observe, that if p = 2 , i.e. A{x, eu) = ueu the system becomes the classical (linear) Stokes 
equation. The non-linear dependence of |ett| is due to the fact, that in the case of non- 
Newtonian fluids the viscosity changes with the velocity (more explicit by the shear rate). 
Power-law or Carreau type fluids are therefore widely used among engineers. For a more 
detailed discussion of the physical model we refer to 1211 Elj and references therein. 

Please observe, that the weak formulation of (HID and (HS]) are the same on dive^ence 
free testfunctions. Indeed, Assumption 11.11 and implies that there exists a uniquqj u G 
LP{{0,T),V,%{n))nL°°{{0,T),L^{n)), such that 

/ udt^ + A{x, Du)D^ dx dt + / u{y,T)^{T,x) — uo{x)^{0,x)dx 
In Jn 

[ j fidx dt for C e {Co“ (H), div^ = 0 }, 

Jo Jn 

for every right hand side / G LP'{[0,T), (FQ^(^v(f^))*) and every uq G L^(H). 

Although our techniques are shaped for the reduced Stokes power law model (jl.Sp they 
are of relevance for the full Navier-Stokes power law model that includes the convective term 

Ut — div(i^|ett -|- fif~‘^eu) + [Vtt]ri -|- Vvr = /. 

It is clear, that for sufficiently large p the convective term can be treated as right hand side. 

For instationary power law fluids little regularity is known. In fact, the results mainly 
relay on classical difference quotients. If they are applied in time direction they imply 
L°°(( 0 ,T),L^(H)) estimates for ut, see [3] (and Proposition 14.111 . In [T7j the authors used 
space and time difference quotient to deduce by embedding Holder continuity for gradients, 
in two space dimensions. See also m where some fractional derivatives of ut are shown in 
the non-degenerate case and in two space dimensions. 

The results present here might very well be an important step to be able to transfer sta¬ 
tionary theory to the instationary case. However, namely in case of three space dimension the 
additional regularity ut G ((0, T), L^(H)) and ut G L^((0, T), A^3:’^(M’*) is of independent 
value. 



Let us briefly discuss the techniques we will use and develop. For the estimate of the 
fractional time derivative we construct backward-forward quotients of the type 

ut{t,x) - f^utit + s,x)ds 

This quotient technique was developed in [laiii] but used on different types of PDF. 


^Here Fo'div(P) := Iff G Co“(P) : divy = 
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For the spatial estimates we will use what we call “Queer Quotients”, which are a mixed 
in time-space difference quotient of of type (ILZI). By doing so we can use the variational 
structure of the elliptic part of the system. This can be combined with the estimate for pure 
in time fractional derivatives already obtained to get fractional differentiability in space for 
ut- Interestingly the estimates are mainly driven by the term ut\ although the ellipticity 
assumptions are needed to treat the term A{-,Du). 

Indeed, the positive terms in our estimates will always be derived from the term ut', not 
like more common from both terms of the right hand side. 


Remark 1.3. AssumMion M. IV (e). which might at first seem the most restrictive assumption 
is actually a consequence of (c). As A{x,-) is differentiable we deduce from (c), that 


1 

¥ 


{A{x, Q + hP) — A{x^ Q)) ■ hP - 


V{Q + hP) - V{Q) 


h 


By letting h ^ 0 we get DqA{x,Q)P iSi P ^ \DqV{Q)\‘^\P\‘^ . Now (e) follows by calculating 
DqV{Q). 


2. Assumptions and Main results 


Our estimates allow to assume growth condition related to a convex Orlicz function. 

Assumption 2.1. Let </? G oo), (0, oo)) 0 oo), [0, oo)), be a convex, ^^'(O) = 0 

and (p'{t) > 0 , for t > 0 function, that holds (p''{t)t‘^ ip{t), for t > 0 . 

This includes (p{t) = for 1 < p < oo and also (p{t) = max{tP,t'^}, for two different 
exponents. Excluded are tlog{t) on the one end and e* on the other end. 

Assumption 2.2. Let (p hold Assumption \2.1[ Additional we assume the following growth 
conditions for x € M” and Q,P £ and // > 0. 

(a) F is measurable and F{x, •) G C'°’^(R"'^^) 0 C'^(R”^'^ \ {0}). 

(b) For A{x,Q) := DqF{x, Q), \ip'{\Q\ + ij)\Q\ < A{x,Q) ■ Q and \ A{x,Q)\ < Klp"{\Q\ + 

h)\QV ^_ 

(c) {A{x,Q) - A{x,P)){Q - P) ~ \V{Q) - E(P)|^ for V{Q) := 

(d) \D^A{x,Q)\ < ap'{\Q\ + fj,). 

(e) PqA(x,Q)P®P~ |P|V"(IQI+/«)■ 

(f) P(x,0)=0. 

In case we consider the symmetric gradient, we need additionally that 

(g) If Q is symmetric, then A{x, Q) is also symmetric. 

The constants used in this assumption are called the characteristics of F. 

In case of ip{t) = F we have the typical p-growth, which is included in the assumption 
above. 

Assumption 12.21 include the following case of Uhlenbeck structure, that F{x,Q) = ip{\Q\). 
Then we have Aj (x, Q) ip'{\Q\)^^. This particular case of Uhlenbeck structure is the sub¬ 
ject of study by an increasing number of mathematicians nowadays. Many physical relevant 
models are of Orlicz growth and not p-growth. For example Carreau type fluids. 

We introduce p*(t) := suPq>o(®^ “ is the conjugate Orlicz function. The following 

inequality of Young type holds (by definition) 

( 2 . 1 ) ah < p*{a) + (p{b), for all a,b £ [ 0 , oo). 

The assumption (p''{f)f^ ~ p{f) implies that its conjugate p* has analogous properties, 
see mm- The space L‘^(0) := {/ : J(^ 99 (|/|) < oo} is a Banach space endowed with the 
Luxembourg norm. Moreover, we define 

Uo'’‘"(L!) and := {5 G Co-(L!) : div <7 = 
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For the definition and properties of and for the analysis on divergence free 

Orlicz spaces we refer to mm- The following estimate is important for us. It can be found 
in uni Lemma 2.4], 

( 2 . 2 ) 

Together with Remark 11.31 this implies that whenever it is well defined 
(2.3) \^^V{Du)\^ ~ + \Q\)\d,Duf. 

Now, standard existence theory implies that if Assumption 12.21 holds, then there exists a 
unique solution of (|1.4p or (|1.5p for uq G L^(n) and / G ([0, T), (Fq^’'^)*). Nevertheless, 
because we need more regularity of the solutions and to keep the paper self contained, we 
include an existence result (Proposition 14.IP . 

We can now state our results. The first theorem concerns fractional derivatives of ut in 
time direction. 


Theorem 2.3. Let F hold Assumption [T^ or Assumvtion \2.‘A Let uq G L^(n) and f G 
L'^* ([0, T), (FQ^’‘^(n))*). Let u be a solution of (|1.4I) or (|1.5p . If additionally f G 1T^’^((0, T),L‘^{Ll)), 
then Ut & — a), L‘^{Ll)) and 

K 

for every a > 0. The constant K only depends on the characteristics of F and the regularity 
off. 


In case of a bounded domain Q, the following estimate is available 

“ ffl {WT* i\f\)\\L^Q t) + \\f\\w^’H{0,T),L^m)- 
The second theorem concerns fractional derivatives of ut in space direction. 

Theorem 2.4. Let F hold Assumvtion \1.1\ or Assumvtion IH.HI on the wholespace M"". Let 
uq G L^(]R”) and f G L‘^*([0,T), (Fq^’‘^(M"’))*). Let u be a solution of ()1.4I) or (jl.Sp . 

If additionally f G ft G IP‘{Qt) and V/ G L‘^* {Qt), then, 

K 

The constant K only depends on the characteristics of F and the regularity of f. 

Theorem 12.31 and Theorem 12.41 and (|1.2p imply the following corollary 
Corollary 2.5. With the same hypothesis as in Theorem \2.4\ we find for u a solution of 
or (fT31) . that ut G AfiofiQr) and therefore ut G L^ (Qt), for q G [2,2-^^), by (fT^ . 

n-\- 2 

In case F has p-growth (i.e. Assumption II.II holds) we can refine the regularity of ut in 
dependence of p. If p > 2 we have to restrict to space periodic solutions. I.e. solutions in 

-7-IIV'llp 

(M"") = b G (^^(IR”) : g{a + x) = f{x) for some o G and / 5 = o| 

This is a closed subspace of F^’^PM”"), therefore a solution of dm) exists. The solution 
then is periodic and it is enough to estimate it over one period, which is a bounded cube. 
These solutions are artificial, however, they are useful to introduce new techniques and new 
observations. In our case we introduce periodic solutions to empathize possible extensions of 
regularity for ut which are to some extend inherited in our techniques. 
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Theorem 2.6. With the same hypothesis as in Theorem \2.4\ but replacing Assumption \2.^ 
by Assumution M. 1\ In case the system (dD is in terms of the full gradient we have: 


(a) p >2 and u E V^er (M”), 

ut E For fl = i + 1} 

Moreover, Ut E T), for q E [l, min , by embedding. 

” 2i3 

(b) l<p<2 andu^V^'P{mA) ut ^ L‘^{{a,T - a),M^ 

.p + 2- {2-p) 


*)), for 


jd = min 


. fl ri 
in I —, max | - 


3 


2p 


6 n . . i 1 

’4 “ 


For m this improvement does not work. However, it is possible to show ut has a /3- 
fractional symmetric gradient. See Remark I5.61 

We conclude the section with the following corollary. It collects the differentiability regu¬ 
larity on Ut and Sobolev embedding o. 


Corollary 2.7. With the same hypothesis as in Theorem. iP.ti] ii>e have Ui € r^^((0, T), L^(M'‘))n 

Kc(QT)r:L]J(0,T),LUn)- If 


(2D) q E [1,4), a E [1,3) for p e (1,4] or q e [1, |],o E [1, forp£ (l,oo) 
g E [1,3), a E [1, |)/or p E [|, 4] or g E [1, ^ 

(nD) q E [l,#^),a E [1,2^^^) for p E \2 - i 


(nD) g E _ 

P E (2 - ^ 


. . ... ..a € [1, /orp E (l,oo), 

p e [2 - 2,4] or q e [l,;^],a £ [1, JJI] fm 


3. Preliminary 

3.1. Notation. We will write p ~ /i if there exist two constants c, C, such that eg < h < Cg. 
For a set A with finite measure, we use 



A 


We define for g : R”" -A R^ 

'■= (g(xi,..,Xi + h,...,Xn) - g(xi,..,Xi,...,Xn)) 

and 

D^iia) ■= j^{g{xi,-,Xi + h, ...,Xn) - g{xi,..,Xi,...,Xn)). 

For an arbitrary direction v E 5”“^ we define 

■■= {g{x + hv) -g{x)) 

and 

F>^{9) ■= ^{gix + hv) - g{x)). 

Moreover, we define the ’’Queer Quotient” as 

:= {9{xi, ..,Xi -h h, ...,Xj + h, ...,Xn) - g(xi, ..,Xi, ...,Xj + h, ...,x„)), 

and 

Dyi{g) := ^{g{xi,--iXi + h, ...,Xj + h, ...,Xn) - g{xi, ..,Xi, ...,Xj + h, ...,Xn)), 
for arbitrary i,j. 

Remark, that in case of u being a solution of m, difference-quotients of u are divergence- 
free and therefore suitable test functions for (HSl). 
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We will use the following partial summation rules, without further reference 


[ A^{f){x)g{x)dx = f fix)A ^{g){x)dx 
Ja J a-\-h 


rb-\-h 


ra-\-h 


+ / fix)g{x -h)- f{x)g{x), 


[ A ^{f){x)g{x)dx = [ f{x)A^{g){x)dx- f f{x)g{x)+ ( f{x + h)g{x) 

Ja Ja Jb—h Ja—h 


and 


D‘ 


pb pb-\-h pa+h 

\f){x)g{x)dx = - f{x)D~'^{g){x)dx + -l- f{x)g{x-h)--l- f{x)g{x). 

J a+h J b J a 


Please notice also the following cancellation property we will use 

(3.1) 

pb pb-\-h pb—h pb-\-h pa-\-h 

/ f{t + h) - f{t -h)dt= / f{t) dt- f{t) dt= f{t) dt- f{t) dt. 

Ja J a+h J a—h Jb—h J a—h 

Finally the following observation is needed (for the sake of completeness). We find for g G 
VF^’^((a, 6)) and for s < /i 

b h b h ^ ^ b h-\- b 

(3.2) f \{D^g){t)\dt < [ -I \g'{t + s)\dsdt = I f \g'it)\dtds< f \g'{t)\dt. 

Ja Ja J J J a+s Ja 


3.2. Preliminary estimates. We start with the following elementary pointwise estimate for 
real functions. It is, however, essential for all our estimates. We use the following estimate 
which is a direct consequence of [121 Lemma 2.7]. If cp"{t)t'^ ~ (p{t) and ao,ai G M"' we have 

(3.3) vp"(|ai| + |ao|) ~ / (p"{\eao + (1 - 6)ai\)d9. 

Jo 

Observe, that although the estimates in Lemma 13.11 and Corollary 13.21 below are purely 
analytic (i.e. not depending on the PDF). We use the notation of the PDF as we will use it. 

Lemma 3.1. Let F hold Assumption or Ig.gl And Du : 11 —There exists 
constants depending only on the characteristics of F, such that for every /i > 0 the following 
identities hold. 

(a) A{Du){t) ■ D'f{Du){t) + ^\A'lV{Du){t)f < D^F{Du){t) < A{Du){t) ■ D’f{Du){t) + 
f\A^V{Du){t)f, 

(b) A{Du){tyDf\Du){t)-i\Af^V{Du){t)f < Df'^F{Du){t) < A{Du){tyDf\Du){t)+ 
y\Ay^v{Du){t)\\ 

(c) 

A{Dum • ^ {Du{t + h)- Du{t - h)) + ^\A>lV{Dumf - ^\Ay^V{Du){t)f 
<^{F{Duit + h))-F{Duit-h))) 

< A{Du){t) • ^{Du{t + h)- Du{t - h)) + ^\A>iV{Du){t)f + ^\Ai^V{Du){tf 
for any fixed value. I.e. ^(•) := A{x, •) and F(-) := F{x, •). 
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Proof. We use the monotonicity of Assumption 12.21 to find 

DiF(-,Du){t) = — [ —F{sDu(t + h) + (1 — s)Du{t)) ds 
h Jq ds 

= f A{sDu{t + h) + {1 — s)Du{t)) ■ Di{Du){t) ds 
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M 


= / {A{sDu{t + h) + {1 — s)Du{t)) — A{Du{t))) ■ D^{Du){t) ds 
Jo 

+ 


f A{Du{t)) ■ D^{Du){t) ds 
Jo 

1 

= / —(AiDuit) + siDuit + h) — Duit))) — AiDuit))) ■ siDuit + h) — Du(t)) ds 

Jo 

+ A{Du{t)) ■ D^{Du){t) 

We estimate using Assumption I2.2lf d) and ()2.2p to gain 


— [A{Du{t) + sA'jDu{t)) — A{Du{t))^ ■ s{A^{Du){t)) ds 


f 


+ \sA’j{Du){t)\ + \Du{t)\)s^\A’jDuf ds. 


Because \sA^{Du){t)\ + \Du{t)\ = s{\A^{Du){t)\ + \Du{t)\) + (1 — s)|L)tt(t)| we use (|3.3I) . 
Assumption 12.21 and ()2.2p to find 


ip'fn + |sAj^(Z)M)(t)| + \Duit)\)s‘^\A^Du\^ ds ~ ^\A'lV{Du)it)\ 


sh 


h' 


This implies the first estimate. 

We redo the argument to hnd that 


D^^F{Du){t) = DfF{Du){t - h) 

= — [ -^F{sDu{t) + {1 — s)Du{t — h)) ds 
h Jq ds 

= f A{sDu{t) + {1 — s)Du{t — h)) ■ D^{Du){t — h) ds 
Jo 

= [ {A{Du{t) + (1 - s){-A^Du{t - h))) - A{Du{t))) ■ D^{Du){t - h) ds 
Jo 

+ f A{Du{t)) ■ D^{Du){t) ds 
Jo 

Jq (1 - - (1 - s)A^Du{t - h)) - A{Du{t))) • (1 - s){-A'f{Du){t - h)) ds 

+ A{Du{t)) ■ D^{Du){t - h) 

We estimate again using Assumption I2.2f dl and ([27 


J + (1 - s)A}Du(t - h)) - A(Du{t))) ■ (1 - s)(-Aj‘(Dti)(t)) ds 

~ jF-^y{A+\{l-i>)Ar’'{Du){t)\ + |Ou(()|)(l - sf\Ai'‘Dufds 
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Because |sA^ ^{Du){t)\ + \Du{t)\ = s(|Aj ^{Du){t)\ + \Du{t)\) + (1 — s)\Du{t)\ we can use 
(j3.3|) . (12.2p and Assumption 12.21 to find 

~ i|A,-V(Du)(()|t 

This implies the second estimate. 

For the third estimate we use 

A{Du{t)) ■ {Du{t + h) — Du{t — h)) 

= A{Du{t)){Du{t + h) — Du{t)) — A{Du{t)) ■ {Du{t — h) — Du{t)), 

therefore 

A{Du{t)) ■ — [Du{t + h) — Du{t — Kj) = - [A{Du{t))D^ {Du{t)) + A{Du{t))D^^{Du{t ))), 

which implies by estimate one and two the third estimate because 

D'iF{Du){t) + D^^F{Du){t) = ^ (^F{Du{t + h) - F{Du{t - h))^ . 

□ 


Corollary 3.2. For every h > 0 and for v G 5” the following identities hold. 

(a) A{Du){z) ■ D^{Du){z) + {Du){z)\ < Dl^F{Du){z) < A{Du){z) • D’^{Du){z) + 

f\A>iV{Du){z)\, 

(b) A{Du){zyDf’^{Du)iz)-f^\A^V{Du)iz)\ < Df’^F{Du){z) < A{Du){zyDf’^{Du)iz)+ 
y\A^;;V{Du){z)\, 

(c) 


A{Du){z) ■ ^{Du{z + hv) - Du{x - hv)) + ^\A^’^V{Du){z)\ - ^\A^V{Du){z)\ 


2h' 


1 


< — [F{Du{z + hv)) — F{Du{z — hv))) 

< A{Du){z) ■ ^{Du{z + hv) - Du{x - hv)) + ^\Ay'^V{Du){z)\ + ^\A’);V{Du){z)\. 


2h' 


The constants only depend on the characteristics of F. 


Proof. As the arguments in the proof of Lemma 13.11 works for arbitrary directions the proof 
is line by line the same. □ 

We will need the following analytic Lemma, which was hinted to us by L. Diening. A 
similar estimate can be found in [51 Lemma 12]. It is an integral characterization of Nikolskij 
spaces. 


Lemma 3.3. Let 1 < q < oo. If f G L'^{[a, b + H]). If 


sup 

H>h>0 


(3.4) 

then f G M^’'^{[a,b]), moreover 


/ \fit + s) - fit)\ds 
Jo 


Li{[a,b+h]) 


< h^K, 


sup |||/(t + /i)-/(i)|||L,([a,b]) <3h^K. 
f>h>0 
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Remark 3.4. Clearly, the reverse direetion is true as well. Simply because 
/ \f{t + s) - f{t)\ds < sup \f{t + s)-f{t)\. 

Jo h>s>0 

Previously one used the slightly weaker but more obvious statement that the bound of 
implied that f € b)), for all ft < a [HI Appendix]. 

Proof of Lemma I,'?. ,4 We estimate 

\f{t + h) — f{t)\ = -/ f{t + h) — f{t + h + s)ds + -l- f {t + h + s) — f{t) ds 
Jo Jo 

^ \fit + h) — f{t + h + s)\ds + -[- \f{t) — f{t + h + s)\ds 
Jo Jo 

rh i‘2h 

\f{t + h) — f{t + h + s)\ds + 2-j- \f{t) — f{t + a)\da. 

Jo Jo 


This implies, for 2h < H 

+ h) - /(•)lllL<J([a,6]) 


< 


/ \f{- + s) - fi-)\ds 
Jo 


+ 2 


([a-\-h,b-\-h]) 


p2h 

t \f{- + s) - f{-)\ds 
Jo 


< 3h^K. 


L<i{[a,b]) 


□ 


4. Estimates in time direction 


We introduce some natural estimates, which are closely connected to the existence theory 
for (II.4p and (|1.5I) . 

Proposition 4.1. Let Assumption\KEhold, 0 < 2a < T and f G L‘^{QT)riL^* 

Then there exists a unique solution u of or (ESI), such that 

f [ \ut\'^ dt dx + sup f F{Du)dx<— f f ip{\Du\) + \f\'^ dx dt. 

Ja Jq (a,T)JQ ^Jo Jo. 

If additionally f G iy^’^((0.r); 

\dtV {Du)\'^ dx dt + sup f \ut\'^ dx < c f ( \dtf\^dxdt + — f f \ut\‘^ dx. 

tG\2a.T] JQ Ja Jq. ^JaJfl 


f 

ha JQ, 


The above estimates are standard. Formally they can be derived by using ut and d^u 
as testfunctions. These are solenoidal, such that the pressure does not interfere. In case of 
p-structure the second estimate can be found in |23l Theorem 8.1]. For the Orlicz setting 
we refer to |3]. As our assumptions are not immediately covered by these references. For 
this reason and to keep the paper self-contained we include a proof of Proposition 14.11 in the 
appendix below. 

To simplify notation we will assume, that u holds (|1.4I) or (II3D on a larger interval say 
{—A, T + A) with initial data at the point —A. We can then assume in the following (by 
using ([32])), that for all h G (0, A\ 

(4.1) f [\D^V {Du)'\‘ + \ut\^ dx dt + sup f \ut\‘^ + F (Du) dx <—, 

Jo Jn te[o,r] Jn A 

Here K depends on / by the right hand side of Proposition l4.ll If LI is bounded one can esti¬ 
mate Jq (p{\Du\) dx dt < (p*{\f\), by testing with u, Poincare and Young’s inequality 
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for Orlicz functions (12.11) . Then 

r / ^*i\f\)+\f\^+\ftfdxdt. 

Jo Jn 

Remark 4.2. Please observe, that in case of p-growth and 1 < p < 2, this implies that 
Dut € LP{{a,T) x fl). As Young’s inequality and ()2.3p imply 

\dtDuf‘ = {\Duf‘-^\dtDu\^)^Du\^^ < \Duf‘-^\dtDu\^ + \Duf ~ \dtV{Du)\ + \Du\p. 


We will need to use ut as a testfunction for our general assumptions. The next lemma 
justihes it. 


Lemma 4.3. Let u be a solution to (naD or (HSl), such that gu holds. Then ut can be 
used as a testfunction to (HH) or (USD, in the sense that 




■ Dut + Uf Ut dx dt 



f • Ut- 


Proof. We start with the following observation which follows from ()2.3p and Assumption 12.21 


r [ \Ai-,Du)\\Dut\ 

o CL kJ 


< < c 

< c 

< c 


f [ (p''{pi + \Du\)\Du\\Dut\dz 
J a JQ. 

f [ ip"{pt + \Du\){\Du\^ + \Dut\^) dz 

o CL kJ 

f f ip{\Du\)+ \dtV{Du)\‘^ dz. 

J a Jo. 


This implies that A{-,Du) ■ Dut G L^{{a,T) x Ll). Next we will approximate ut. We use 
the approximation D^u, for 0 < h < a which is a testfunction to the system on the interval 
[a,T — a\. We then have pointwisely, that 


A{x, Du{t, x)) ■ D^{Du){t, x) —>■ A{x, Du{t, x)) ■ Dut{t, x) for —>• 0 


and almost every (t, x). To be able to use Lebesgues convergence theorem we have to provide 
a majorant. Let t G [a,T — o]. Please observe, that by ()2.2p we have that 

\D^V {Du){t)\^ ~ ip"{p + max {I Hu (t) I, \Du{t + h)\})\D^{Du){t)'^. 


If \Du{t + h)\ < \Du{t)\, we estimate using the previous and Young’s inequality 
lA(Du(t)) ■ Hf(Hu)| < ip"(p + |Hu(t)|)|Hu||Hf(Hu)(t)| < cip{\Du\) + c\D^V{Du){t)\^. 

In case \Du{t + h)\ > \Du{t)\, we estimate using Assumption 12.21 and the last estimate to get 
\A{Du{t)) ■ D^{Du)\ < \ {A{Du{t)) - A{Du{t + h))) ■ D^{Du)\ + \A{Du{t + h)) ■ D'1 {Du)\ 

< c\D’lV{Du){t)\^ + c\D^V{Du){t + h)\^ + cp{\Du{t + h)\). 


All together we have the majorant c{\D^V{Du){t)f + \D^V{Du)(t + h)f + ip{\Du(t + h)\) + 
p{\Du{t)\)) which is converging in L^((a,T — a) x H). As we assume / G L'^{Qt) and as we 
know by ()4.1I1 that ut G L'^{Qt), we find that 



A{-,Du) ■ DD^u + Ut ■ D^u dx dt 



f • Dfu, 


converges to the right limit equation with /i —)■ 0. 


□ 


The following Proposition is the main effort to prove Theorem 12.31 
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Proposition 4.4. Let the assumptions of Proposition \4-l\ be satisfied, such that holds. 
Then ut G Moreover, for every A > h > 0 we have 


h rT-A-h 


fl 


0 JA+h 


iAfK)(t)r 

h 


, , , , K 
dx dtdsda < 

~ A 


where K depends (linearly) on the right hand side of (j4.1l) and on the characteristics of F. 

Proof. Let A > h > 0 and a G [0,h]. We use the testfunction ^fQA'j~^Af{ut){t)ds on 
[a + h,T — a — h]. Let us briefly show we are allowed to use this test function for every 
positive h. We calculate 


/ 


^A){ut){t) ds = 2ut — u{t + h) — u{t — h), 


the functions u{t + h),u{t — h) £ V^’‘^{Qt) and therefore admissible. The function ut is 
admissible by Lemma 14.31 

Therefore we may apply the testfunction to our system. Partial summation implies 


(/) + (//) + (///):=!£ 

1 


1 

^h 



T—a—h 

\Af{ut){t)\'^ dxdtds 
0 JQ Ja-\-h 

a+/i 

utAf{ut) dt dx — 

0 J Q, J a-\-h—s 
T—a—h 



T—a—h 

utA){ut) dtdxds 

Q. JT—a—h—s 



—a—n pn 

A{Du){t)D-j- 2ut{t) — ut{t + s) — ut{t — s) ds dx dt 
h Jo 


= -/ 

hh 


J a+h J 0 

h r pT—a—h 

At{f)At{ut)dt 
0 J a-^h 

T—a—h 


pa-\-h pi—a—n 

+ fA){ut)dt— / f Al{ut) dtdx ds =-. {IV) 

J a-\-h—s JT—a—h—s 


We start by estimating {II) with Young’s inequality 


ph pa+h p pi —a—n 

|(L/)| <-f + / \ut{t)\‘^ + \ut{t + ,s)\'^ dt +-I \ut{t)\'^ + \ut{t + s)f dt dx ds 

Jo J a Jo, J T—2h—a 


T—a—h 


< c sup / 

{a,2a+h)U{T—2h—a,T—a) JQ. 


\ut\^ dx. 


Next we estimate {IV) by Young’s inequality and similar as before 


\{IV)\<5{I) + c,\^j^ 


T—a—h 


\Al{f)\dxdtds + c sup 

0 JnJa+h (a,2a+h)U{T—2/i—a,T—a) 


[ kl" + l/r 

Jn 


dx. 


We divide {HI) into 


= -H 


pT-a-h 

A{Du){t)Dut{t) dx dt 

IQ. J a-\-h 

pT-a-h 


^ r pJ—n—n 

A{Du{t)){Du{t + h) — Du{t — h)) ds dx dt = {III)i + {III)-^ 
h Jq J a+h 


urn = y 


pT — a—h 2 P 

dtF{Du{t)) dx dt =/ F{Du{T — a — h) — F{Du{a + h)) dx 

fQJ a+h ^ Jo. 
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To estimate (1 11)2 we need to use Lemma l3.1 


{ 111)2 = 4 / / {F{Du{t - h)) - F{Du{t + h))) 

^ Jn Ja+h 

- ^\AT>^V{Du){t)\ + ^\A’^V{Du){t)fdxdt 


> —c 

+ 


II 


T—a—h 


a+h 
(i-\-‘2,h 


A^V{Du){t) 


h 


dx dt 


rT-a 


[ [ F{Du{s)) dx ds — I I F{Du{s))dxds 

h ^Ja Jn JT-a-2h Jq 


= —c 


II 


T—a—h 


Q J a+h 


IS>iV{Du){t) 


h 


2 , ra+2h 


[ F{Du(s)) 
Jn 


dx dt 

rT-a 


dxds-J- / F{Du{s))dxds 

T—a—2h J n 


where we used m- This implies 


(I) <c\{II)\+c{\{IV)\-6{I)) + c 


T—a—h 


J a+h 


+V{Du){t) 


h 


dx dt 


(4.2) 


2 ra+2h p 

+ TT / F{Du{s)) dx ds — F{Du{a + h)) dx ds 
"-J a Jn 

+ —-/ f F{Du{T — a — h)) — F{Du{s)) dx ds. 

hj T-a-2h Jn 


The terms in consideration have the factor which looks to much. However, by an inte¬ 
gration with respect to a, there is an additional cancellation effect which gives the “correct” 
order in h. This procedure is an important moment in the proof. 

We integrate a over the interval [0, A] to find 


(4.3) 


pA pa+2h 

h Jo T a 


/ F{Du{s)) dx ds — F{Du{a + h)) dx ds da 
'n 

2 p pA p2h pA+h 

= — / / J F{Du{a + s)) dads — / F{Du{a))dadx 

f^JnJo Jo Jh 

p p2h pA+s pA+h 

= — / -/ / F{Du{a)) da ds — / F{Du{a))dadx 

“ JnJo Js Jh 

p p2h ph pA+h 

~ J J I F{Du{a))da— j F{Du{a)) dadsdx 


'0 Js 

r2h 


lA+s 


f -f F{Du{a)) da — -f F{Du{a)) dal ds dx 

JnJo y s J A+s ^ 


< sup / F{Du) dx. 

{0,h)\j{A+h,A+2h) J 
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The second difference of F is estimated analogous, such that 

r*A I'T—a p 


f 


T-a-2h Jn 


F{Du{T — a — h)) — F{Du{s)) dx ds da 

T-h 


(4.4) 


-I p pr—n pzn p 

T - F{Du{a)) da + f / 

Jn Jr-A-h Jo Jt 


1 f 

^ JnJo 


pT'—A—2h-\-s 

iT-A-h 


2h pT—2h-\-s 

F{Du{a)) da ds dx 

T-A-2h+s 

T—2h+s 

F{Du{a)) da ds dx 

T-h 


F{Du{a))da+ f f 

Jn Jt- 


< sup / 

(T-A-2h,T-A-h)U(T-h,T) Jn 

Combining (|4.2li with (14.811 and (14.4p gives 

rT-A-h . 


F{Du) dx. 


fl 


ff I r 

Jo Jo Jn J a+h 


0 JA+h 


m 


h 


< 


+ 


sup / 

(0,A+2/i)U(r-A-2/i,T) Jn 

h pT-h r 

2 


dxdtds < 
F{Du) dx + 


T—a—h 




h 


dx dt ds da 


sup 

(0,2h+A)U(T-2Ii-A,T) 


[ Wtf 

Jn 


dx 


-if 

^Jo Jh 

'// 


Jn 

T—h A h 


A.l{f)\ dx dt ds + sup / 

(o,2/x+A)u(r-2h-A,r) Jn 

2 


utV + \ f\dx 


A>lV{Du)it) 


h 


dxdt] < c 


K 

A 


by Sobolev embedding and (|4.1I) . 


□ 


Proof of Theorem \2.,J[ We imply ProDosition l4.4l on the interval on [|, T] for A = |. Lemma[331 
then concludes the proof. □ 


5. Estimates in space direction 


In this section we derive estimates for mixed derivatives (in time and space) for the whole 
space situation. Therefore If = M” and Qt = (OjT) x M"’. 

The aim of this section is to show fractional derivatives in all directions including slanted 
ones in space-time. Naturally we will assume more regularity of the data in space direction. 
The following proposition is the space-analogue of Proposition 14.11 It can be derived by 
formally testing with —Au 

Proposition 5.1. Let Assumvtion \2.S\ hold. 0 < 2a < T and f G ((0, T), (I/^’‘^(n))*). If 
additionally V/ € {Qt) we have for u the solution of ()1.4p or m- 

[ f \W {Du)\'^ dx dt + sup f \Vu\‘^dx<c f f (p*{\V f\) + (f{\Vu\) dx dt. 

J a te[a,r] Jr^ Jo 4r" 

The proof of the above proposition for (II.4p follows by the difference quotient technique. 
Basically the technique of [51 Theorem 6.1] can be applied to the parabolic situation with 
general F. For the convenience of the reader we include a proof in the appendix. 

As before to simplify notation we shift the time interval by A, such that we can assume in 
the following 


u 


+ sup 
tG[0,T] 


/ 

Jv 


\D^V{Du)\ +\dA{Du)\ +\uA + ip{\Du\)dxdt 


\Vu\^ + \uA + F{Du)dx < 


(5.1) 
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for all 0 < /i < A. Where K depends on / by Proposition 15.11 and ()4.1I1 . The next step is to 
redo the argument before in space directions. 

Lemma 5.2. Let u be a solution to or dLSl), and let the assumptions on f be such that 
dSII) holds. Then 

rh rT 


[ [ ^xi{ut){x)Ki{uxi){x)dxds 

I^Jo Jo 


K 


Where ^ is the right hand side of m- 

Proof. We use the test function A~f .{uxf){x) ds. This implies for almost every t G 

[o,r] 

(/) + {II) + {III) := U [ A^,^{ut){x)A^,^{ux,){x)dxds 
i^Jo Jw^ 

+ J I A{Du){x) ■ D-f 2uxi{x) - Uxi{x + sci) - Uxi{x - sci) ds{x) dx 
h Jr™ Jo 

= [ Al^{f)-A%^{ux,)dxds=-.{IV) 

'XJo 4r™ 

The term {IV) can be estimated 

\{IV)\<f [ g,*{\D^,J{x)\) + g,{\Du\)dx. 

Jo Jr'* 

The terms {II) can be estimated exactly as in the proof of Proposition 14.41 

{II) = — [ A{Du{x)) ■ di{Du{x))2 +f A{Du{x)) ■ {Du{x) — Du{x + hci)) dx 
“ Jr™ u Jr™ 

9 r 

A{Du{x)) ■ {Du{x — hci) — Du{x)) dx 


+ 


2 

h 


f diF{Du{x)) dx — — f A{Du{x)){D^.{Du{x)) + D^I^{Du{x)) dx. 
Jr™ h JjRn 


The first term vanishes, the second can be estimated by Corollary 13.21 
|(//)| < / {DlF{Du{x)) + Df^F{Du){x)) dx 


+ c 


I 

Jv 


A-I^V{Du) 


h 


< c 


K 

A' 


□ 


We can combine Proposition 14.41 and the last lemma to get diagonal directions. 


Lemma 5.3. Let u be a solution to (d or (d, and let the assumption on f be such that 
(d holds. Then for every h,A>0 


pA ph p 

Jo Jo J a+h 


T—a-\-h 


V^{ut) ■ Ay{d/^u) 

h 


dx dt ds da 


K 

< — 
- A 


where ^ depends (linearly) on the right hand side of (j5.1|) by the characteristics of F. 
Proof. We will use the testfunction 
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which is admissible by the same argument as in the proof of Proposition 14.41 We integrate 
over [a + h,T — a — h] x M"" and find 

(I) + {II) + {III) 


-I 


pT—a—h 

/ {ut{t,x))A~^Ay<{d/^u)dxdtds 
Ja-\-h 

pT—a—h 

/ A{Du){t, x) ■ D{d/^u){t,x) dx dt 
I a-\-h 

pT—a—h 


—a—n pn 

A{Du){t, x)D+ dy^u{t + s, X + ejs) + dy^u{t — s,x — Cis) ds dx dt 
h Jo 

1 r pT—a—h 

— j- / / f A~JA\{d/yu) dsdxdt =-. {IV) 

hJo Jmx Ja+h ^ 


We start by estimating (!) with partial summation 




T —a—h—s 

{Ay,ut){t, x)Ay<{d^u) dx dt ds 


' a+h 


\ r po,-\-h 

+ TT / / ut{t + s,x + eis)Ay^{d^u){t,x) dtdx 

^Jo Jr.'^ Ja-\-h—s 

= {I)l + {1)2 - {1)3. 


pT—a—h 

/ utA^y^ddt dx ds 

IT—a—h—s 


We find that 


We estimate 


(/)l + {II) + {III) = -{1)2 + {1)3 + {IV) 


|-(/)2 + (/)3| < C sup 

tG[a,a-\-2h]\j[T—a—2h,T—a] 


[ Inti’ll-\Du\" 
JR^ 


dx, 


which can be estimated by ()5.1I) (after taking the supremum of some terms over time). 

We estimate {IV) by 

1 C / pT—a—h—s 

^Jo Jr^ \Ja-\-h 

pa-\-h pT—a—h \ 

+ / f{t + s,x + eis)Ay{ut + Uxi){t,x)dt- f{t,x)Ay{ut + Uxi){t,x)dt\dxds, 
J a-\-h—s JT—a—h—s / 

where the terms depending on u can again be estimated by (|5.ip . This implies that 

p ph pT—a—h—s 

\{iv)\< f / \Dy.f\‘^ds + \Du\‘^+ \ut\^dxdt 
-/Jo Ja-\-h 


(5.2) 


+ 2 
+ 2 

+ 2 


// 

If 


a-\-2h p pa+2h 

\ut\^ + \Du\^ dx dt + 2 / -h Ifl'^dxdt 


I, J 


a+h 

T—a—h 


T—a—2h 
h pT—a—h—s 


p pi —a—n 

ut\^ + \Du\^ dx dt + 2 / -j- \f\‘^dxdt 

Jr^J T-a-2h 

pn pi—a—h—s 

j- / \Dyf\‘^ds + \Du\'^ + \ut\'^dxdt 
Jo J a-\-h 

sup f \ut\‘^+ \Duf+ \ff dx 

[T—a—2h,T—a) Jr‘^ 


(a,a+2/i)U(T—a—2 
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The terms {II) and [HI) will be estimated as in the analogous estimates before. 

rT—a—h 


or —a—ri 

{II) = — / / d/'F{Du)dxdt 

h Jbx Ja+h 


The term {HI) we need to use Corollary 13.2 

rT—a—h 


F{Du{T — a — h)) — F{Du{a + h)) dx 


\ f —a—Ai 

(///) = —— 5 - / / A{Du){t,x)D{u{t + h,x + eih) — u{t — h,x — eih)) ds dx dt 

n Ja+h 


2 r pT—a—h 

— T 2 / {F{Du{t — h,x — eih)) — F{Du{t + h,x + eih))) 
“ JR" Ja+h 


- ^\A^W{Du){t) \ + j^\AyV{Du){t)\^ dxdt 


n 

J a JW 


> —c 


2 ra+2h r n pi 

+ Tf / F{Du{s))dxds —I F{Du{s)) dx ds. 
^J a J’SX f^J T-a-2h iR" 

Therefore, we integrate a over (0, A) and have 


AyV{Du 


^ • h 

2 


dx dt 


pT-a 


(■A pA pT-a P A’^.V{Du) 

(/)ida|</ l-(/)2 + (/)3| + |(/C)|+c/ / ^ 

3 Jo Ja JR" " 


dx dt da 


+ 


+ 


l*A 2 pci-\-2h 

Jo ^J a 

pA r\ pT—a 

i If 


^J T-a- 


F{Du{s)) — F{Du{a + h)) dx ds da 
I F{Du{s)) — F{Du{T — a — h)) dx ds da 


2h . 


By (14.31) and (14.41) we gain 


/ (/)i da\ < c sup 

lo lg[0,A+2/i]U[r-A-2h,r] 

r*/i pT—h—s 


f + \Du\'^dx 
jR" 

P ptl pi —ti — s 

F j- \Dy^f\'^ds + \Duf + \ut\'^dxdsdt 

Jr^Jo Jh 

sup f \ut\‘^ + \Du\'^ + \f\'^dx 

A+2h)U{T-A-2h,T) JR" 

pT P a^.V{Du) ^ r 

/ / —- dxdt+ sup / (p{\Du\) 

Jo JR" ^ (0,A+2/i)U(r-A-2/i,T) JO 


dx dt + sup 

(o,A+2/i)u(r-A-2/i,: 

All terms can either be estimated by (|5.1|) such that we gain the result. 
Theorem Ea follows from the following proposition and Lemma 13.31 


dx. 


□ 


Proposition 5.4. Let F hold Assumvtion \l.l\ or As,sumvtion \2.l^ on the wholespace M”". Let 
uq G L^(R”) and f G L‘^*([0,T), (Pq^’‘^(M"’))*). Let u be a solution of ()1.4I) or (jl.Sp . 

If additionally f G bT^’^((0, T), n) and V/ G L^*{Qt), then, 

JO JA 


m JA+h JO 


^fh 


A 
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for every A > h > 0. The eonstant K only depends on the eharacteristies of F and the 
regularity assumptions of f. 

Proof. We start with a few basic calculations. We reformulate 

= ut{t, X + eis) — ut{t, x) = ut{t + s, X + Cjs) - ut{t, x) + ut{t, x + set) - ut{t + s, x + scj) 
= Ay{ut){t,x) + Afut{x + sci). 

Therefore, we find 


\Al.ut\ < \Ay{ut){t,x)\ + \Afut{x + sei)\. 

The last term can eventually be estimated by Proposition l4.4[ Therefore, once we can estimate 
the ’’Queer Quotient”, we gain the result. He is estimated by 

|A^(uQ(t,x)p = Ay{ut){t,x) • Ay{ut + Uxi){t,x) - A^(uQ(t,x) • Ayu^iit^x) = Sq- Si 

The term Sq can be estimated by Lemma 15.31 
We analyse Si 

Si = Ay{ut){t,x) ■ Ay{u^.) + Ay{ut){t,x) ■ Afiu^.) 

(5.3) + Af(uQ(t,x) • Al.{u^.) +At{ut){t,x) ■ Afin^J 

!= All T AI 2 T AT^ T 71^4 


This implies by Proposition 14.41 and Lemma 15.31 

r-h rT-A-h p 

/ I A®, (nt)(t)|^ dx dt ds 


pn pi—/ 

Jo J A+h 


ffj 
/77 

Jo Jo Ja 


< 


< 


pT—a—h 


' a+h 
pA ph pT-A-h 


+ 

+ 


A®, (ui) (t) I ^ dx dt ds da 
A®(ttQ(t)|^ dxdtds 

A+h JK.^ 
pA ph pT—a—h 

TT Ay{ut){t) ■ Ay{d/-u) dx dt ds da 

7 0 ■/ 0 J a+h 


/7 

*do do 


hK 


A ph pT—a—h 

All + AI 2 + M 3 + M 4 dx dt ds da 

0 J a+h 

A ph pT—a—h ph pT—h 

Ml + M 2 dx dt ds da 

a+h 


pA ph pi 

UO Jo 7aH 


+ 


ph pT—h p 

j- / |Mi + M 2 I dx dt ds. 

do Jh dR" 


Ml can be estimated by Lemma 15.21 


(5.4) 


pA ph pi 

Jo Jo J a- 


f'A ph pT—a—h 

All dx dt ds da 
1—h 


< 


Kh 

~A~' 


We are left to estimate the integrals with M 2 , M 3 and M 4 . Unlike the other terms they 
can only be estimated by less, i.e. by +hK/A. We estimate the integral of M 3 with Holder 
and Young. 


f / -f \ ■ At,.{u+\dsdtdx 

w h dR" do do 

r A\At{u+t,x+ 

Jr^ Jo 


<csup / |idrt| dx + c 
[o,r] dR’^ 


Jo 


ds dt dx. 
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The integral with M 4 can be estimated analogous: 

rT-h ph 


f f -f \Af{ut){t,x) ■ Af{uxi)\dsdtdx 
VhjR^Jo Jo 

f in -L f f'^\^t{ut){t,x)\^ 

ap / \Du\ dx + c ' ' ' 

,T] JR" 


< c sup 
[ 0 / 


/ If 

JR"JO JO 


ds dt dx. 


To estimate M 2 we use partial summation and proceed as before 
(5.5) 

1 

\/hJo Jo J 


< 


pT—a—h 

/ A^^.{ut){t,x) ■ Af{uxi) dt ds dx da 
Ja-\-h 

2 r rh pT—a—h 

-^= 4 - I A^^{ut){t,x) ■ A~^Uxidtdsdxda 
vhJo JR" JO Ja+h+s 


+ 


A rh 


VhJo Jo 


pT—a—s /*a+/i+s 

/ A^^.{ut){t,x) ■ {uxi)dt - Al.{ut){t,x) ■ {uxjjdtdsdxda 

/ T —a—h J a+/i 

=:(/)+ {II) 

The estimate on {II) can be estimated optimally 

(5.6) {II)<2'/h sup / iDu]"^ + \utf dx. 

[0,2h]n[T-2h,T] JR" 

The term (/) is again estimated by 

rT-h /■/»|A|(^,)(t,x)|2 




(/) < c sup f I Du\ ^ dx + c f f J- 
[0,T] JR" JR" JO JO 


h 


■ ds dt dx. 


□ 


Next we will prove Theorem 12.61 In the following we assume power-law structure. I.e. 
that Assumption 11.11 hold for (11.41) which is assumed to be in terms of full gradients (i.e. 
Du = Vn). 


Proof of Theorem \2.6[ . We will use the following Sobolev embedding 


r f 2(n + l) 

/ / \V{Du)\ "-1 dxdt< 

Jo JR" 

We will also use the following estimate in case n > 3 


\VV{Du)\'^+ \dtV{Du)fdxdt] < 00 . 


n+1 
n —1 


r [ |Vn| 

JO JR" 


^ dxdt < 


< 


TT 2 

[ [ (\Vu\P^dx') " [ f|Vu| 

JO JR" V / JR" V 

7T 2 

/ ( \V{Du)\^^ dx) dt sup 

JO VjR" / te(o ,7 

/ [ \'^V{Du)\‘^ dx dt sup 

JO JR" te(o,T 


_4 n \ n 

"2 dx dt 


|Vuf dx 


|Vuf dx 


Where we used Sobolev embedding. (By Korn’s inequality the same estimates hold also 
in case of symmetric gradients). Observe, that the above estimate implies Vu G L^{Qt)^ 
where q := max {; P + „; ^• Ifl <P< best integrability of Vu will 

be 2, and Theorem 12.41 can not be improved. In case p > min 2 — 1 we enter the proof of 

ProDosition l5. 41 below (15.411 . As mentioned there we only need to estimate the terms involving 
M 2 , M 3 , M 4 (defined in (|5.3p . as all other terms can be estimated by hK/A. 
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We start by estimating M 4 . We define 2/3 = min{^ + 0, 1}, where 9 will be fixed at the 
end of the proof. Young’s inequality implies. 

Y| A|{u,) ■ A;(a,.)| < A| Af (k,)|G Af(u„)l" 

The first term can be estimated by Proposition 14.41 The second term is a mixed derivative, 
(a) Case p G (max{l, 2 — ^}, 2]. By Young’s inequality and the definition of Vi{Du) 

0-2 I 


1 




V(|Vu(t+ s)| + |Vrt(t)|)P-2 /i2 


(2-p)0 


< (|Vu(t)| + |Vu(t + s)|)^(|Vw(t + s)| + |Vu(t)|)^^ + c\DlViiVu)it)\^ = (I) + {II) 


The function {II) is integrable. By Young’s inequality we find 

2-ep 2-ep 

{!) < c|Vu(t)| 1-9 + c|Vrt(t + s)| 1-9 
which is integrable if = <?• Calculations imply, that we can choose 


9 = max 


j-pT 2 - (2 -p)n 

I ^ 



Observe, that 0 > ^ for all p G [max{1,2 — ^},2]. As /3 = 
that /3 = ^ for p G [2—^,2]. 

(b) Case p G [2, 00 ) we fix 0 = | and estimate 


min{^, I + |}, we have 




e 


which is locally integrable. If u is space periodic, so is ut and for functions of this 
type local integrable functions are also globally integrable. This implies that /3 = ^ 
for pG [2,4] if u G Yp^ef(M"). 

The term [Maj can be estimated in the very same way. On the term M 2 we have to apply the 
partial summation as is done in ()5.5I) : there the term {II) can be estimated optimally (see 
(|5.6p and (/) analogous to the term M 4 above. □ 


Remark 5.5. We will give some explanation why large p decrease the differentiability of ut- 
Let us consider the model case dni). If p > 2 , the quantity V{Du) is differentiable and not 
Du itself. By (|2.3p one observes that 

J \dtV{Du)\'^ ~ J \Duf-^\dtDuf < 00 . 

Therefore, dtDu is only integrable with respect to the weight Theorem, \2 . hi sa,us that 

/3 —I for p 00 , which would be Proposition |5..^| again. We conjecture, that for general 

1 

degenerate F when p 00 , the best possible order of mixed derivative is ^f^ut. 


Remark 5.6. Let us discuss the case of symmetric gradients in the framework of Theo- 
rem 1.2.61 Only at the estimates of Mi,..., M 4 we needed the restriction. This is due to the 
fact that VY(Vrt) does not necessarily exist. Only W{eu) is an L'^-function. Therefore, by 
the proof above it is only possible to get analog bounds on 

KM) + ^ 

hh 

in the framework of symmetric gradients. 
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6. Appendix 


Proof of Proposition \4-l\ We briefly recall that and (|1.5p have a unique solution with 
additional regularity. 

A simple and standard way to prove this is based on Ritz-Galerkin Approximations. 

m 

Um{t,x) = 

i=i 

with smooth Ansatz functions E Wq''^{PI) nL^(n); which are linear independent and 

whose linear hull is dense in Wq’‘^( 0) n LP‘{Vl). We choose 

m 

x^ = {v£ n L\n)\v = 

and E X^, such that tt™ —>■ uq strongly in L‘^{Q). The approximate equation reads 

+ (^(•, = {f,Pk), Um{0) = U™ 


for all 1 < /c < m and is equivalent to a system of ODE for the cj*. 

In the existence interval [0, r) we may use the test function u'^ and obtain a uniform 
estimate 


f f \Umfdxdt+ f F{Dum)dx\l < f f \f\'^dxdt 
Jo Jn Jn Jo Jn 

uniformly in m E N. Hence the solution can be extended to the full interval [0, T], Moreover, 
the following estimate holds for a. e. r E (a, T]. 


rp T ^ 

(6.1) / / \u'^fdxdt+ sup / F{Dum){T)x < / / \f'^dxdt + j- / F{Dum)dxdt. 

Ja J D. rE[a,T]./r2 Jo Jq, J Jq, 


Due to convexity the Ritz-Galerkin approximation is equivalent to the minimum problem 


( 6 . 2 ) 



v)dxdt+ f f F{Dum)dxdt< 

Jo Jn 



I 

Jn 


F{Dv) + f{u. 


v) dx dt, 


for all V E Xm and all r E [0, T]. 

By weak compactness we find a subsequence Um u and u'^ u' in and Vttm ^ Xu 
in L'^. In (|6.2p we perform an integration with respect to r E [ti,ti + h] and h small, then 
(IS2I) can be rewritten as 


ti+h 


oil f I u'^v dx dt + [ f F{Dum)dxdtdT 

^ J Jn Jo Jn Jo Jn 


ti+h 


< 


-/ / / F(Dv) + f {um — v) dx dtdr + f dx. 

J Jo Jn Jn 


Now we may pass to the limit m —>■ oo and due to lower semi continuity in the weak topology 
we find 


(6.3) 


u) dx dtdr 


\u\^{T)dx— f f u'vdxdt+ f f F{Du) 

Jo Jn Jo Jn 

f f F{Dv) + f{um — v) dx dtdr + f |uo|^dx. 

Jo Jn Jn 
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for a.e. r and v G {i}) n-L°°((0, T), L^(Q). This can be reformulated into 


/7”‘< 


u — v)dxdt + 



0 Jo, 


F{Du) dx dt < 



0 Jn 


F{Dv) + f{um — v) dx dt 


and a weak solution is constructed. Additionally, (jh.ip holds in the limit and implies the first 
estimate in Proposition 14.11 Ibv Assumption 12.2p . 

We proceed by taking a sequence of Vm G such that Vvm —^ Vu strongly. By taking 
Vm as a testfunction in (16.2h and passing to the limit we find. 


lim sup / / 

m Jo Jn 

This implies (by strict convexity) 


F{Dum) dx dtdr < 


>0 Jn 


F{Du) dx dtdr. 


lim 


/ F(Dum) dx dtdr = / / F(Du) dx dtdr and therefore Dum ^ Du stiongly. 

Jn Jo Jn 


Furthermore, differentiating the Ritz Galerkin equations with respect to t and testing with 
u'^ we obtain by (e) of Assumption 12.21 

1 ^ , f f \d^v{Vu^)\‘^dxdt< 


(6.4) 


dx\l + 


\ftf + \u'^f dxdt. 


'a J 


fa J Q. 


We gain a subsequence dtV{Dum) dtV. However, by the strong convergence of Dum 
we find V = V{Du). Now passing to the limit in (16.4p implies the desired estimates after 
integration over a. □ 


Proof of Provosition 15. it For the sake of simplicity we take F{x,Q) = F{Q) omitting the 
lower order terms. By Assumption [221(c) we have that 

D^,^{A{Du)) ■ D^,^{Du) ~ \D^,y{Du)\\ 

We take the testfunction —D^jD^Ju)), this implies by partial summation and Young’s in¬ 
equality (12.ip 

/ »)-(|Di./l) + A|Of(«)l)- 

JR" 

By passing with h —0 we get the desired estimate. □ 
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